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Abstract; We study the Stokes problem in a bounded planar domain 17 with a friction type boundary condition that 
switches between a slip and no-slip stage. Unlike our previous work El, in the present paper the threshold value may 
depend on the velocity field. Besides the usual velocity-pressure formulation, we introduce an alternative formulation with 
three Lagrange multipliers which allows a more flexible treatment of the impermeability condition as well as optimum 
design problems with cost functions depending on the shear and/or normal stress. Our main goal is to determine under 
which conditions concerning smoothness of 17, solutions to the Stokes system depend continuously on variations of 17. 
Having this result at our disposal, we easily prove the existence of a solution to optimal shape design problems for a large 
class of cost functionals. 
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1 Introduction 

This paper analyses one property of the Stokes system dehned in 17 C with a slip type boundary condition, namely the 
continuous dependence of its solutions on the shape of 17. This property plays the crucial role in the existence analysis of 
optimal shape design problems. The no-slip boundary condition, i.e. the vanishing velocity on the boundary, is widely used 
in practice. It characterizes the adhesion of the fluid on the solid wall. This condition is acceptable for small velocities and 
on a macroscopic level. On the other hand, there are many situations (flow of the fluid on hydrophobic surfaces, polymer 
melts flow, problems with multiple interfaces, micro/nanofluidics etc.) where the slip of the fluid occurs. To get a more 
realistic model, the slip has to be taken into account. For the physical justification of different types of slip conditions we 
refer to Qo) and 121 . The mathematical analysis of the Stokes and Navier-Stokes system with the slip and leak boundary 
conditions has been done in 121 and extended to non-stationary problems in IS). The regularity of solutions to the Stokes 
system with slip and leak boundary conditions has been established in ifTTIl . In HI the stick-slip condition is considered as 
an implicit constitutive equation on the boundary and the existence of weak solutions to Bingham and Navier-Stokes fluids 
is proven. 

Shape optimization involving fluid models with slip boundary conditions as the state problem is of a great practical 
importance. Slip boundary conditions affect the velocity profile and hence the velocity gradient of the fluid in the vicinity 
of the wall. The velocity gradient is an important factor in the transformation of the mechanical energy to heat, the 
process representing the energy loss. Shape optimization of the interior of hydraulic elements may reduce the velocity 
gradient resulting in energy savings. In ||6l a class of shape optimization problems for the Stokes system with the threshold 
boundary conditions involving a priori given slip bound has been studied. The existence result for the continuous setting 
of the problem and convergence analysis for appropriate discretizations of the continuous model have been established. 

Nevertheless it is known from experiments that the slip bound may depend on the solution itself, e.g. on values of the 
tangential component of the velocity. The aim of this paper is to extend the existing stability results to this type of the slip 
boundary condition. Besides the standard velocity-pressure formulation used in ||6l we present a new weak formulation 
adding another two Lagrange multipliers; one releasing the impermeability condition and the other regularizing the non¬ 
smooth slip functional. This new formulation turns out to be useful in numerical solution of this problem. Moreover, 
it enables us to approximate directly the normal and shear stress and to use these quantities as arguments of appropriate 
objective functionals to control the stress distribution along the slip part of the boundary. 

The paper is organized as follows: in Section|2]we present the velocity-pressure formulation of the Stokes system with 
a solution dependent slip bound. Using fixed point arguments we prove that such problem has at least one solution for 
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any slip bound represented by a continuous, positive and bounded from above function g. If in addition, g is Lipschitz 
continuous with modulus of Lipschitz continuity sufficiently small, then the solution is unique. Section deals with a 
four-field formulation of the problem whose solution is represented by the velocity u, pressure p, normal, tangential shear 
stress a'^, and a'^, respectively. In Section |4] we prove that the graph of the respective generally multi-valued solution 
mappings considered as a function of the shape of the slip part of the boundary, is closed in an appropriate topology. On 
the basis of these results the existence of solutions to a class of optimal shape design problems will be proven in Section|5] 

2 The velocity-pressure formulation of the problem 

Unlike in, where the slip bound was given, the present paper deals with a more general case, namely the slip bound will 
be a function of the tangential velocity. 

Let O C be a bounded domain with the Lipschitz boundary dfl. The slip boundary conditions are prescribed on an 
open, non-empty part S of the boundary and the no-slip condition on L = \ S', L ^0; 


Alt + Vp = f 

in U, 

(2.1a) 

div It = 0 

in U, 

(2.1b) 

It = 0 

on F, 

(2.1c) 

Itjy = 0 

on S, 

(2. Id) 

kk<5(KI) 

on S, 

(2.1e) 

Ur ^ 0 ^ \a'^\= g( Mt ) & 3A > 0 : Ut = —Xa'^ 

on S. 

(2.1f) 


Here u = (mi, W 2 ) is the velocity field, p is the pressure and f is the external force. Further, zz, r denote the unit outward 
normal, and tangential vector to SU, respectively. If a S is a vector then := a ■ v>, Ur '■= a ■ t is its normal, and the 
tangential component on dil, respectively. Finally, cr’’ := (f^)^ stands for the shear stress and g : K+ — > R.+ is a given 
slip bound function. By a classical solution of (12. Il l we mean any couple of sufficiently smooth functions {u,p) satisfying 
the differential equations and the boundary conditions in (12.11) . 

To give the weak formulation of (12. Il l we shall need the following function sets: 


V(fl) = {r; S (iF^(U))^| r; = 0 on F, = 0 on S'}, (2.2) 

V^div(^^) = {tt G V(r2)| divD = 0 a.e. in U}, (2.3) 

Lg(U) = {geL2(U)| [ q = 0}, (2.4) 

Jn 

L+(S) = {v? e L^{S)\ p>0 a.e. on S}, (2.5) 

= {ip £ T^(S)| 3z; G z; = 0onr: v = p on S}, (2.6) 

= {p£ <p>0 a.e. on Sj. (2.7) 


Remark 1. Ifv € V (U) and S £ then it is readily seen that z;.r|s £ 

From now on we shall suppose that S £ C^'^. 

The trace space is equipped with the norm 

ll^lli/2 s = inf 

V-y—ip 

where w[ip) £ Y (U) is the solution to 

Aw{ip) = 0 in U, ' 
w{ip) = 0 on F, 
w„{<p) = 0 on S, 

Wr{p) = Lp on S.^ 

Further we introduce the following forms: 

a{u,v)= / Vm : Vt), b{v,q) = / gdivv, j{p,Vr) = / g{p)\vr\, 

Jn Jn Js 

u,v £ q £ p £ h]^^{S). ( 2 . 8 ) 
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We shall assume that g : K.+ —)> R+ is continuous and there exist positive constants gmin < 9max such that 

9min ^ !?(■) — 9max tit ]R_|_. 

The weak formulation of (12, 11 1 reads as follows; 

Find {u,p) £ VX Lq{^) such that 

Mv £V (n) : a{u, v — u) — b{v — u,p) 

> 

+j{\uT\,Vr) -j{\Ur\,Ur) >{f,V- u)o,n, 


(2.9) 


(V) 


Vg G : b{u,q) = 0. 


We will show that under the above mentioned assumptions on g, has at least one solution for any f £ (L^(f}))^. 
To this end we use the weak variant of Schauder’s hxed point theorem jSl. 

1 /2 

For a given function p £ (S) we consider the auxiliary problem: 

Find {u'^,p'^) £ V(fl) x Lglfl) such that 


Vv £V{fl) : a{u’^,v — u'^) — b{v — u'^,p‘^) 

> 

Vr) - j{>f, O >{f,V- u‘^)o,Q, 


(vn 


Vg G : b{u'^, q) = 0. 


1 /2 

We know that for every p £ Flj^ 

l|V«^|lo,o + lb^llo.a<c(ll/ll 


(S') there exists a unique solution {u^ 

o,n + II5 MIIoo,s) ^ c, 


ofdzgi and it satishes (see 0): 


( 2 . 10 ) 


where c := c(/, gmax) is a positive constant. 


Let us dehne the mapping ik : h]^'^{S) — F[]^'^{S) by 
T'((p) = |u^| on S. 

Then (|^ is equivalent to the problem of hnding a hxed point of ik in F[]^‘^{S). 


Theorem 1. The mapping ik has the following properties: 

(i) ik(i?) C B, where B = {p £ h]^‘^{S)\ ||<p ||]^/2 S — andc is the constant from (12.101 1. 

1 /2 

(ii) Ik is weakly continuous in FI^ (S), i.e. 

Pk^p in pk,p £ hI^^S) ^ in 

Proof The property (i) follows immediately from 

IIKIIIi/2.s<IKIIi/2.s<IIv«‘"IIo,o<c, 

making use of (12.10b . 

Let {u^,p^) denote the solution to pk ^ p vet iT^/^(S). Since the sequence {(m^,p^)} is bounded in V(fl) x 

Lq^^I), there exists a subsequence (denoted by the index k') such that 

^ It in p^ ^pin Ll{n), fc'—(■ oo. 

It is easy to show that {u,p) is the solution of (17^‘^b . i.e. {u,p) = (it‘^,p‘^). Indeed, 

limsup a{u^ ,v — ) < a{u, v — u), 

k' —¥oo 

{f,v-u''')o,n^{f,v-u)o^Q, * (2.11) 

b(v — jP’^ ) ^ b{v — u,p) \/v£V{n).^ 
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It remains to show that 


f 9iv’k')i\vr\- \u^ \) ^ [ giip)i\vr\-\ur\), k'^ oo. (2.12) 

Js Js 

To prove (I2.121 i. we use that 

ifk' ^ </3in => (fik" —>■ ^ a.e. on S' g{(pk") gi^) a.e. on S, (2.13) 

where {(fk"} denotes a subsequence of {ipk'}- Moreover, 

^ K| in 2.2(5-). (2.14) 

Clearly (12.13b and (12.14b imply (12.12b . From (12.1 lb and (12.12b it follows that {u,p) is a solution to Since this 

solution is unique, then 

^ (m,p) weakly in (22^(fl))2 x Ll(yt,), fc —>• oo, 
i.e. (m,p) = {u‘^,p'^). Finally 

^ in ^ ^ |m^| in (i2^(f2))2 ^ |m^| ^ \u ^\in H^/‘^{S) 

'^{p) in i2^/2(5) 


proving)//). □ 

From the weak variant of Schauder’s theorem and Theorem[T]it follows that there exists at least one fixed point of T* in 
i2y^(S) and thus at least one solution to (TPb . 

Remark 2. Problem (17^‘^b is well posed also for p G L^(S). Similarly as in the previous theorem, one can show that 
T' : L^(S) L'^{S) is continuous: 

Pk pinL^{S), pk,p G T+(S) => '^(pk) '^(p) inL^{S). 

Next we shall study under which conditions, problem (|3 has a unique solution. 

Theorem 2. In addition to ( 12.9b . let g : R+ —^ R+ be Lipschitz continuous in K+.' 

3L > 0 : \g{xi) - g{x 2 )\ < L\xi - X 2 \ 'ixi,X 2 G M+. (2.15) 

Then T* : (S) —>■ L^(S) is Lipschitz continuous, as well: 

||T'((^i) - ^'(v52)|Io,s ^ \\pi - (/52|Io,s G L\^S), (2.16) 

where L is from (12.15b and c is the norm of the trace mapping tr : V (Ll) —>■ L‘^{S), trv = Vr, assuming that V (Lt) is 
equipped with the norm \ ■ |i.n- 

Proof Let pi,p 2 G L'^{S) and {u\p'^) be solutions to / = 1, 2. Then it® G Vdiv(i2) and 

a{u\v - u') + j(K|,'!;r) - j(KI>0 >{f,v- M*)o,n Vv G Vdw{il). 

By a standard technique we obtain: 

- «^) < / ig{pi) - g{^ 2 )){\ul\ - 1 ) 

Js 

<|05) L \\pi - V?2|lo,S ||l«rl - l4l||o.S ^ \\Pl - ‘F2|lo,S l^r “ Wr||o.S 
< cL \\pi - p2\\o,s 1“^ - w^li.n- 

From this and the fact that 


||T'((/?i) - T'(v?2)|Io,s = lll^rl - |w?|||o.5 <c\u^-u^\i,n 

we obtain (12.16b . 


□ 
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Corollary 1 . If L < 1/c^ then 'I' is a contraction in L'^{S). Its unique fixed point belongs to the set Bi = {p ^ 
^+('5')! \\V’\\o,s ^ cc}, since 

IllWrIllo.S = \\Ut\\o,S ^ < CC, 

as follows from (12.101 . Consequently, the method of successive approximations 

•Co G L\{S) given, 

Ck+i ■= '^{•Pk), fc = 0,1,..., 

converges in L^{S) to the unique fixed point z of ^ in 

Ck ^ z in LF‘{S), '^{z) = zonS. (2.18) 

In fact, z € 11^^“^ (S). Indeed, let u^, u be the first component of the solution to and {V^), respectively. It is 

sufficient to show that z = |wr | € From (12.171 and (12.181 it follows: 

\u^ - M|i,n < cL \\tpk - -2:|lo,s ^0, oo, 
so that 

Wck - |Mr||lo,S = llkrl “ l^-rlllo.S - 11^^ “ '“'^llo.S - /c -5> OO. 

From this and (12.181 it follows that z = \ut \ on S. It is also readily seen that 
^ p in L^lyS) 

and {u,p) is the solution to (V). 


3 Four-field formulation of (\V‘^h and (TPl) 


The pressure p in the velocity-pressure formulation introduced in the previous section is the Lagrange multiplier associated 
with the incompressibility condition in ft. This section presents another formulation involving two additional Lagrange 
multipliers a'^, a'^ defined on S releasing the impermeability condition = 0 on S, and regularizing the non-differentiable 
functional j. To this end we shall need the additional function spaces: 


W{n) = {veH^{n)\v = 0onT}, (3.1) 

w{ii) = w{n) X w{n), (3.2) 

77-1/2(5') = (771/2(5))' (dual of 771 / 2 ( 5 ))^ (3 3) 

7fi/2(5) = 77^/2(5) X 771/2(5'), (3.4) 

77-1/2(5) = (771/2(5))'. (3.5) 


If/r = ( 71 , 72 ) G 77-1/2(5), if = (( 71 ,^ 2 ) G 771/2(5) then 

{f,c) '■= {Ci,‘Pi) + ( 72 ,‘ 72 )- 
Since 5 G Ci’i, the mapping 

tr : r; !->■ (v„,Vr), where v„ = t ;|5 ■ i/, Vr = ^is • r, 
maps W(n} onto 771 / 2 (5). ^ e then 

{n,tTv) := + {p",Vr). 

Analogously to the previous section the space is equipped with the norm 

llv’lli /2 s = inf = l«i(¥’)|i,n, (3.6) 

tr v—(p 
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where w{ip) in (13.6b solves: 

Aw{(p) = 0 in O, 
w{(p) = 0 on r, 
trw{(p) = ip on S. 

The standard dual norm in is given by 

M-1/2.S = sup = sup 

t>GW(0) 


<pP0 


trvpO 


(M.tru) 


where 


ll/2,S 


is defined by (13.6b . One can introduce another norm on H {S), namely 
(/r,trt!) 


IIMII-i/ 2,S= sup 

vew(n) l^’li.n 

v^O 

It is known IH that 

M- 1 / 2 .S = IImILi/ 2 .s V/r e 

To regularize the functional j{ip, •) we introduce the bounded convex set K(p) defined by 

K{p) = {/i- e L^{S)\ ImI < 9{p) a.e. on 5}, G hI^^S). 

It is readily seen that 

j(p,Vr)= / g{p)\Vr\= sup / Vr ■ 

Js iJ.^GK(v)Js 

Hence 


j{<P,VT) > {p"',Vr)o,S V/r^GiT((^). 


(3.7) 


(3.8) 


(3.9) 


T\\& four-field formulation of p G reads as follows: 

Find{u,p,a'',a^) G W{n) X Ll{n) X x K{p) s.t. 

Vu G VI^(O) : a{u,v) - b[v,p) - {a'',vf) - (cr"^,■u^)o,s = (/,i^)o,n, 

VgGLg(O): 6 (m, g) = 0, 

G 77-1/2(^) : = 

^K{p)\ + cr^,Wr)o.s < 0. 




Suppose that (IA7‘^b has a solution. In what follows we give its interpretation. Clearly from (IA7‘^b o 3 we see that u G 
^div(i7), where Vdiv(^7) is defined by (12.3b . Using test functions v gV (U) in (lAd’^b j we get 


a{u,v- u) - b{v - u,p) - {a^,Vr - Ur)o,S = {f,v- M)o.n Vu G V(U). 

From it follows that 

-{cr'" ,Ur)o,S = sup {p^,Ur)o,S = j{<P,Ur), 

which together with (13.9b yields 

-(cr^,Ur -Ur)o,S <j{.9^,Vr) -j(p,Ur). 

From this and (13.10b we obtain: 

a{u, V - u) - b{v - u,p) + j{p, Vr) - j{p, Ur) >{f,V- M)o,n Vu G y (U), 
i.e. the couple {u,p) G V (U) x Tq(U) solves iV^i . The formal application of Green’s formula to (IA7‘^b ] gives: 


(3.10) 


cr‘" = —p + 


du 

dv 


and cr’^ = 


du 

du 


on S. 
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Theorem 3. Problem (1^4*^!) has a unique solution (u,p, a'^,a'^) for any ip S In addition, the couple {u,p) 

solves (17^‘^b . 

Proof. To prove the existence and uniqueness of a solution to it is sufficient to show that the bilinear form c{v, q, pf := 
—b{v, q) — (/r, trt;) satisfies the LBB-condition. It is well-known (see e.g. ||4l) that 


37 >0: sup > 7lkllo,a W e Lg(n) 

v^O 


and from (I3.8l l we have 


(/r,trt!) 

sup - 

vew(n) |u|i,a 
vpo 




Then there exists a constant 7 > 0 such that 

> 7(ll9llo.o + M- 1 / 2 .S) V(g,/r) e Lg(fl) x 

vew{Q) Fli,n 

vpo 

as follows from Theorem 3.1 in □ 


(M) 


Any quadruplet (it, p, ,a'^) is said to be a solution of the problem with the solution dependent coefficient of friction 
if it solves dSfJ with (fi = |mt-| on S: 

{u,p,a'',a^) G W{n) X Lg(fl) x x K{\ur\) s.t. 

Mv G TF(n) : a{u,v) - b{v,p) - {a'' ,vf) - (cr^,Vr)o,S = (f,v)o,n, 

VqeLlifl): b{u,q) = 0, 

yp’'G : (/i",M,) = 0, 

G Ki\ur\) : {p^ + a^,Ur)o,S<0- 

On the basis of the results of Section|2]and Theorem|3]we arrive at the following theorem. 

Theorem 4. Problem iA4i has a solution. In addition, if g satisfies (12.151 1 with L as in Corollary\I\ then the solution is 
unique. 


4 Stability of solutions with respect to boundary variations 

The aim of this section is to show that solutions to and (lAfb depend continuously on the shape of O. We shall suppose 
that only the part S of dVl where the slip conditions are prescribed, is subject to variations. In addition, for the sake of 
simplicity of our presentation we shall assume that S is represented by the graph of a function a which belongs to an 
appropriate class Uad- Here and in what follows Uad will be defined by 

Uad = {aG C'P^([0,1])| 0 < amin <a< amax in [0,1], < Cj, j = 1,2 a.e. in (0,1)}, (4.1) 

where amin, amax and Cj > 0, j = 1, 2 are given. With any a G Uad we associate the domain 

H(a) = {ixi,X 2 ) G M^l xi G (0,1), X 2 G (a(xi),a;)}, 

where w > 0 is a constant which does not depend on a G Uad- The family of admissible domains consists of all r2(a) with 
a G Uad- We shall also assume that / G (L^(R^))^. 
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4.1 Stability of© 

Let a G Uad be given and denote by {u{a),p{a)) G V (^^(q!)) x LQ{il{a)) a (not necessarily unique) solution to {V{a)) 
defined in 17(a): 


Vi; G V(17(a;)) : aa{u{a),v — u{a)) — ba{v — u{a)),p{a)) ' 

+ jai\Uria)\,Vr) - ja ( |Mt (ct) |, Ur («) ) 

>{f,v- M(a))o,n(a), 

\fq G Ll{n{a)) : ba{u{a),q) = 0, 


{V{a)) 


where Ua, ba and are defined by (12.81 1 on 17 := 17(a), S := S{a). 


Remark 3. Let us note that the condition on L ensuring the uniqueness of the solution to {V{a)) can be chosen to be 
independent of a G Uad- Indeed, the constant c in (12,161 1 can be bounded from above uniformly with respect to a G Uad 
(see Lemma 2.19 in m- 

Let 17 D 17(a) Va G Uad be a hold-all domain and tt^ G C{V (17(a)), (i?Q (17))^) an extension mapping from 17(a) to 
17. Since all 17(a), a G Uad satisfy the uniform cone property, there exists tTq, whose norm can be estimated independently 
of a G Uad- Finally, the upper index stands for the zero extension of functions from 17(a) to 17. 


Theorem 5. There exists a constant c := c{f, g-max) > 0 independent of a G Uad such that 


IKaM(a)|li^o + ||F°(«)|lo.n ^ 


(4.2) 


holds for any solution {u[a),p{aj) to {V{a)). 

Proof The estimate of the first term in (I4.21 i follows from (3.4) in 16). Similarly, the estimate of the second term in (14.21) 
is the same as (3.5) in jh) using that g < gmax in K+. □ 

Let 


Gv ■= {(a, M(a),p(a))| a G Uad, {u{a),p{a)) solves (F’(a))} 
be the graph of the generally multivalued solution mapping 4) : a i—;> {u{a),p(^a)), a G Uad- 
Theorem 6. The graph Q-p is closed in the following sense: 
a„ a in (^^([0,1]), a„,a G Uad, 

(7ra„Mn,P°) ^ («,p) in {Hl[n)f X Ll{n), 

where (a„,M„,p„) := (a„, M(a„),p(a„)) G Qv. 
and hence (a, M|n(a),P|n(a)) G Gv- 


(“|n(a),P|a(a)) solves {V{a)) 


Proof Let v G V (17(a)) be given. From Lemma 3 in ||6l we know that there exist: a sequence {vk}, Vk G {H^ (^))^ and 
a filter of indices {uk} such that 


Vk ^ V in (iF^(17))^, fc —>■ cx) 


(4.3) 


and 


Ufc|n(a„^) G ^^(^(ctnfc))- 

Therefore Vk\Q(^ani^) can be used as a test function in {V{anf))'- 

O-oin^, i'O'nk J Ufc — tl-nk ) ~ (Ufe — linj, , Prik ) 

+ ja„f. (lUnfcrl j Ufer) ~ (|UnfcT |, Unkr) 

— if ,'Ok )o,n(o!„j _), ^ 

Vg G Lo(^(an J) : (M»fc, 9) = 0. 


(4.4) 
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Denote (M(a),p(a)) := (M|n(a))P|n(a))- The fact that it (a) G Vdiv(f^(a)) and the following limit passages can be 
proven exactly as in m- 

limsupaa„ {u„^,Vk - Un^) < aa{u{a), v - u{a)), 

k—>oo 

lim {Vk - = ba{v - u{a),p{a)), 

k—^oo ” 

= {f,v- M(a))o.n(a„j- 

It remains to pass to the limit in the slip terms given by j. 

From (I4.31 t and weak convergence i'^nk) ^ ^ in it follows (see Lemma 2.21, ||5|): 


“nfc|S(a„J 

0 M|S(q) 

o a 1 

in(T^(0,l))" 

^fc|S(a„,) 

O an^: '*^|S(c<) 

o a \ 

and also 



Un^T O Qf„j, 

■= (t^rifclSian^) 


UkT ° Clnj, ■ — (ttfe|S(a„j,) 

■T“"»=)°anfc (r’is(a) -t"; 


o a = Ur o a 
o a = Vr o a 


in L\0,1), 


(4.5) 


using that o an^. =1 t“ o a (uniformly) in [0,1], where stands for the unit tangential vector to 5'(/3), /3 G lAad (see 
ll6l). Consequently, there exists a subsequence of {un^r ° o-n ^} (denoted by the same symbol) such that 


UrikT o ctrik Ur O a a.e. in (0,1). 


(4.6) 


Hence 


(It^rafcT I) ffer) — J O Orifc \')VkT O Onj, dxi 

—^ / 9{\Ur O Q;|)t^T O Q;\/1 + {a')'^dxi = JadWrl, Vt), k —>■ CXD, 

Jo 

making use of (l4.5b T. (14.61 1 and convergence Q!„ —> a in C^([0,1]). The limit passage for the second slip term in (I4.41 l can 
be done in the same way. □ 


4.2 Stability oflTMl) 

Let Uad be defined again by (lO . We keep notation of Subsection l4.ll i.e. H(a), S{a), ja, Ua, be t“ have the same 
meaning. In addition, (, )„ will denote the duality pairing between and iJ^/^(S'(a)). If G H~^/^{S{a)) 

and V G (Tf^(H(a)))^ then 

(/^ , Vi>)a ■— {p ; • ly )a Ver G l^ad 

and similarly for {p'^, Vr)a, P^ G H~^/'^{S{a)). Problem (lAfb formulated on H := n{a) will be denoted by {M{a)). 

Let 

Qm = {(a,M(a),p(a),cr''(a),CT'^(a))| a G Uad, {u{a),p{a),a''{a),a'^{a)) solves {M{a))} 
be the graph of the respective solution mapping. 

Theorem 7. The graph Qm is closed in the following sense: Let 

a„ -)• a in 1]), a„, a G U^d (4.7) 

and 

(7ra„Mn,P°) ^ («,p) in (Tfo(H))^ x Lo(H). (4.8) 

Then also 

^ (cr‘'(a),u^)a, (4.9) 

^ (o-^(a),'yr)o,S(a) (4.10) 

holds for every v G (Tfg where lUn,Pn, CF'^^ <Jh) is a solution to (AI(a„)), n = 1,... In addition, the quadruplet 

(w|n(a),P|n(a)>0'''(«):Cr^(«)) solves (A4(a)). 
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To prove this theorem we shall need the following auxiliary result. 

Lemma 1. Let g,h G L°“((0,1)), g > 0 be such that \h\ < g a.e. in (0,1) and {gn}, gn G L°°{{0, 1)) be a sequence of 
nonnegative functions, gn ^ g o..e. in (0,1). Then there exists a sequence {hn\, |h„| < gn a.e. in (0,1) such that hn ^ h 
a.e. in (0,1). 

Proof. The interval (0,1) will be decomposed as follows: (0,1) = where 

h = {x G (0,1)1 h{x) = 0}, 
h = {xG (0,1)1 h{x) = g{x), g{x) ^ 0}, 
l 3 = {x G (0,1)1 h{x) = -gix), g{x) f 0}, 
h = {xG (0,1)1 h{x) G (0,5(x))}, 

J5 = {a; G (0,1)1 h{x) G {-g{x), 0)}. 

The sequence {h„} dehned by 


0 

on/i, 

9 n 

on/2, 

~ 9 n 

on/3, 

min(/i,p„) 

on /4, 

_max(/i,5„) 

on 4 


has the required properties. □ 


Proof of Theorem^ In virtue of Theoremj^we have to prove (14.91 1 and (14.101) only. We use the formulation of (A^(a„)): 
{Un,Pn,Crn^^n) ^ W {n{an)) X Lg(fl(a„)) X (S (an)) X K{\UnT\S{a„)\) ■ 


yv G (H^oin)) 


2 . 


{o'n,Vn)an =««„(««, It) “ (it, Pn) 

'*^t)o,S(q:„) ~ (/;'f)o,0(a„) I 


Wq G L‘^{n{an)) ■■ &„„(««,<?) = 0, 

V/r'" G i/“i/2(5'(a„)) : {g'', Unu)a„ = 0, 

V/r G iT(|Unr|S(a„) I) • (p T O'ni ^™'r)o,S(a„) 0 , 


(4.11) 


where 


K{\unr\s{a„)\) = iP"' ^ L'^(S{an))\ |/r^| < 5(|M„|S(a„) ''t""!) a.e. in ^(an)}. 

Observe that one can use test functions v G {Hq (O))^ in (14.Ill) , since W{n{f3)) = (Hq (ll))^f2(^) V/3 G Wad- 

Denote u{a) := u\n{a) andp(a) := P|q(„). Letting n —>■ cx) in (14.111) , we obtain: 

lim {(Jn,v„)a,, = aa{u{a),v) - ba{v,p{a)) - (/,t>)o,n(a) - lim «,'yT)o,S(a„)- (4.12) 

n—¥oo n—foo 

Further 


lim (cr 

n—foo 


T 

n’ 


^r)o,S(a„) 


= lim 


/ ■ T' 

JS{a„) 


= lim [ dn(vis(a„} ■ T°''‘) o an^/T+Ja^ dxi, (4.13) 

where = < o a„. Clearly |(T^| < 5(|tW|s(a„)I) a.e. in (0,1) with tW|s(a„) = (Mn|S(a„) • t“") o a„. Hence there 
exists a subsequence of {&„} (denoted by the same symbol) and an element 9'^ G Lf{{Q, 1)) such that 


r in L^((0,1)). 


(4.14) 
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It is easy to show that \9'^\ < 5(ttT(a)|s(Q,)) in (0; 1) making use of (I4.61 l and (14.141) . Since a„ a in 1]) and 

t"” o a„ t“ o a in [0,1], we finally obtain; 


lim «,U^)o,S(a„) 

n—>-cxD ^ 


[ 0^iv\S(a) ■ t“) O a^/l + {a')^ dxi = (O'" ,Vr)o,S{a), 

Jo 


(4.15) 


where 9'^ is defined by 9'^ = 9'^ o a. Thus 9'^ € iT(|Mr(Q!)| 5 (Q) |). 

Since u{a) G Vdiv(^^(Q;)) we have 

ba{u{a),q) = 0\/q € Ll{n{a)), (4.16) 

and 

ifi'', uM))o. = 0 G H-^/\S{a)). (4.17) 


To verify the last inequality (I4.11L . let G iT(|uT-(a)|s(a)l) be arbitrary. Then \ji^\ < (7(|wT(Q;)|g(^)|) a.e. in (0,1). 
From Lemma[T]and the Lebesgue theorem we know that there exists a sequence |/)^| < <7(|tW|s(a„) I) in (0,1) 
such that 

An A^ in T^((0,1)). 

Consequently, 

0 > (/in + finr)o.S(a„) = / (An + K){'^n\S(a^) ’ t“") O \/l + «)2 dxi 

Jo 

f + 9^){u{a)\s(c.) ^r^)oa^/l+~(^ dxi 
Jo 

So far we have shown that there exists a function 9'^ G K{\ur{ct)\s(a)\) such that 
lim (cr((,r’i.)a„ = aa{u{a),v) - ba{v,p{a)) 

n—^oo 

- (('^,WT)o.S(a) - (/,ii)o.n(a) Vt> G (i^/((^7))^ 

(/i^ +6»^,Un(a))o.S(a) <0 V/i^ G iT(|Un(Q;)|S(a)|), 

+ (|4~T^ and (HTTtI i. 

To finish the proof we have to show that 9'^ = a'^(a). Indeed, then (14.91 ) follows from the fact that the right hand side of 
(I4.18b i is equal to (a"(a), Vu)a and (14.101 ) from (14.151 ). We use again Lemma 3 from ||6l: for any v G V(f2(Q;)) there 
exists a sequence {vk}, Vk G and a filter of indices {n^} such that 

-)> Uin (i7o(fl))^ (4.19) 

and 

■Wfc|n(a„j G y(0(anj), (4.20) 

where tJ|o(a) = v. The definition of {A4{an^)) yields; 

b “ (^nfc ) ij^nk j "^fc) ~ i^k^Puk ) ~ (^n^ 5 ffer )o,S(ct„j,) ~ (/) )o,n(a„|^ ) ■ 

Letting k ^ oo and using (14.181 ) we obtain as before; 

0 = aa{u{a),v) - ba{v,p{a)) - (0^, Wr)o,S(a) - (/,'w)o,n(a) = (CT^(a),'yr)o.S(a) “ , Vr)o,S{c,) 

holds for any v €V (fl(a)). Hence <j'^{a) = 0'^. □ 


= id" + S'" ,Ur{a))o^s{a)- 


(4.18) 
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5 Application of the stability property in optimal shape design problems 


On the basis of the results of Section |4] it is easy to prove the existence of solutions to a class of optimal shape design 
problems for systems governed by the Stokes equation with a solution dependent slip bound. 

Let J-p : 073 —R and Jm '■ Qm —R be cost functionals defined on the graphs of the solution mappings corre¬ 
sponding to {V{a)) and a G Uadi respectively. Shape optimization problems with {V{a)), as the state 

relation read as follows: 


Find z* G Gv such that 


> 

Jv{z*)<Jv{z)^zGgv,^ 


(P 73 ) 


where 2 * = {a*,u{a*),p{a*)), z = (a, u{a),p{a)), and 
Find z* G Gm such that '1 

Jm{z*) < Jm{z) 'iz gGm,] 


(Pa^) 


where 2 * = {a*,u{a*),p{a*),a''{a*)), z = (a, tt(a),p(a), cr'^(a), cr"^(a)), respectively. 

To prove the existence of solutions to (|F^ and (|F^ we use compactness and lower semicontinuity arguments. 
Convergence in Gv and Gm will be introduced using the results of Theorem|6]and|2] We say that Zn —> 2 :, Zn, z G Gv 
if (14.7b and (14.8b hold. Analogously, Zn -G z, Zn, z G Gm if (I4.7b - (l4.10b hold. 

From the definition of Uad and the Arzela-Ascoli theorem we see that Uad is a compact subset of . This, together 
with (14.2b and Theorem|6]and|7]proves the following result. 

Theorem 8. The sets Gv <^nd Gm compact with respect to convergences introduced above. 

We say that Jp and JTat are lower semicontinuous functionals on Gv, and 07^1, respectively if 


Zn -S' 2:, Zn, z GGv ^ livCLini Jp{Zn) > Jv{z), 

n—¥oo 


(5.1) 


Zn Z, Zn, Z gGm ^ liminf Jm (zn) > Jm(z). (5.2) 

n—¥oo 

Theorem 9. Let the functionals Jp, Jm satisfy dSB, and (15:2b . respectively. Then there exists a solution to i jP-pb and 



Proof is straightforward. 


Conclusions 

In the first part of this paper we analyzed the mathematical model of the Stokes system with a threshold slip boundary 
condition whose slip bound depends on the solution itself. We used two weak formulations of this problem: the standard 
velocity-pressure formulation and the extended formulation in terms of the velocity, pressure, shear and normal stress. 
We proved the existence of a solution for a large class of functions representing the slip bound and studied under which 
conditions the solution is unique. In the second part of the paper we analyzed how solutions to both weak formulations 
depend on the geometry of the problem, in particular on the shape of the slip part S of the boundary. Using an appropriate 
parametrization of S we proved that the graphs of the respective solution mappings are compact in an appropriate topology. 
This result plays the key role in shape optimization. 
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